Ultrafast sequential charge transfer in a double quantum dot 



O 

(N 

< 

m 



^ ■ 



> 

in 

o 
o 



X 



A. Putaja and E. RasanerQ 
Nanoscience Center, Department of Physics, University of Jyvaskyla, FI-4OOI4 Jyvaskyla, Finland 

(Dated: September 1, 2010) 

We use optimal control theory to construct external electric fields which coherently transfer the 
electronic charge in a double quantum-dot system. Without truncation of the eigenstates we operate 
on desired superpositions of the states in order to prepare the system to a localized state and to 
coherently transfer the charge from one well to another. Within a fixed time interval, the optimal 
processes are shown to occur through several excited states. The obtained yields are generally 
between 99 % and 99.99 % depending on the field constraints, and they are not dramatically affected 
by strict frequency filters which make the fields (e.g., laser pulses) closer to experimental realism. 
Finally we demonstrate that our scheme provides simple access to hundreds of sequential processes 
in charge localization while preserving the high fidelity. 

PACS numbers: 78.67.Hc, 73.21. La, 78.20.Bh, 03.67.Bg 



I. INTRODUCTION 

During the past few years, coherent control of charge in 
double quantum dots (DQDs) has been a subject of active 
experimentali"— and theoretical^^— research. Here one of 
the long-term aims is the design of a solid-state quantum 
computing scheme.-™ It is still to be seen whether the 
optimal control mechanism DQDs turns out to operate 
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through magnetic fields ; ' gate voltages^ or optimized 
laser pulses,'''^ 

Dynamical control of charge in DQDs has been a pop- 
ular application for few-level scheme3^ d°d'^~ '22i (modeling 
DQDs as two-, three-, or four-level systems), which have 
demonstrated ultrafast high-fidelity processes. However, 
a physical DQD has, in principle, infinitely many lev- 
els, and in fast processes a considerable number of states 
might have practical relevance. For example, a two-level 
approximation is exact only in the limit of using an in- 
finitely long resonant continuous wave with an infinitely 
small amplitude. A linear field (bias) is an appealing and 
simple alternative to control charge in DQDs,^ but it is 
not applicable to fast processes. 

With quantum optimal control theor y^°i^^ (OCT) it is 
possible to find optimized external fields driving the sys- 
tem - having an arbitrary number of states - from the 
initial state to the desired target state without any ap- 
proximations, apart from a possible model potential to 
describe the physical apparatus. OCT has been used 
to analyze the general controllability criteria of two- 
dimensional single-electron DQDs and to optimize in- 
terdot charge transfer ii Optimal control of two-electron 
DQDs has been obtained in an extensive work of Nepstad 
at al.^ addressing various control schemes^ and hyper- 
fine interactionsi^ 

In this work we apply OCT to construct external elec- 
tric fields that lead to fast sequential charge transfer pro- 
cesses in single-electron DQDs. To obtain high fidelity we 
operate on the superpositions of the lowest states corre- 
sponding to the charge localization in left or right well. 
We show that hundreds of sequential charge transfer pro- 
cesses can be achieved without a significant loss of the 




FIG. 1: (Color online) Model potential for the quantum dot 
(black solid line), the ground state (green line), the first ex- 
cited state (black dash-dotted line), and their superpositions 
corresponding to left (red dotted line) and right (blue dashed 
line) states. 



yield. To make the experimental production of the ob- 
tained fields more realistic, we cut off the high-frequency 
components already during the optimization procedure. 
The use of such filters does not dramatically affect the 
fidelity. 



II. MODEL 

We use a one-dimensional (ID) model describing a 
single-electron semiconductor DQD. The external poten- 
tial has a form 



{(-f)^(-fr} 



(1) 



in effective atomic units (a.u.), see below. Here d — 6 
is the interdot distance and = 0.5 is the confine- 
ment strength. The potential is visualized in Fig. [1] 
We consider typical GaAs material parameters within 
the effective-mass approximation, i.e., m* = 0.067 and 
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e = 12.7, 

scale as 



El = 



Now, the energies, lengths, and times 
{m* /mo)l{e/eofEh « llmeV, = 



{e/eo)/ {m* /mo)ao ~ 10 nm, and — h/E^ 



60 fs. 



spectively. We emphasize that below the abbreviation 
a.u. refers to these effective atomic units. 

It should be noted that a harmonic potential in Eq. ([T]) 
is, in its two-dimensional (2D) form, a general model 
for realistic semiconductor quantum-dot structuresj^ 
Since the first Coulomb-blockade experiments it has been 
shown that the harmonic model is essentially valid up to 
dozens of electrons confined in the dot, and thus up to 
a large number of levels. The validity is clear, e.g., in 
recent works combining experiments and theory in the 
spin-blockade regime . ^^'^^ The precise energy-level spec- 
trum in a given device can be explicitly obtained through 
single-electron transport experiments, and this informa- 
tion can be utilized to reconstruct the particular form of 
the extenal potential. For example, in Ref. [l^ it was ex- 
plicitly shown that measured energy-level spectrum can 
be well reproduced by a harmonic model potential upon 
slight refinements. Hence, when necessary, Eq. ([T|) can 
be tuned to match a particular device. Regarding the 
results below, the ID model does not yield a qualitative 
difference from a more realistic 2D potential, but it sig- 
nificantly speeds up the calculations. 

Electronic states localized to left and right dots can 
be expressed as superpositions of the two lowest (gerade 
and ungerade) states |0) and |1) as follows: 



\L) = -^(|0)+|1)) 

\R) = ^(|0)- |i)) 



(2) 
(3) 



If the system is prepared in either of the superpositions, 
the occupation probabilities of \L) and \R) oscillate with 
the resonance frequency wqi = i?i — i?o ~ 0.0135 (see 
Ref. Issh . For instance, if the system is first prepared 
at \L), it reaches the state \R) a.t t = T/2 = tt/llIqi ~ 
232.87. As discussed in detail below, we aim at control- 
ling this charge-transfer procedure in an arbitrary way. 



where the (static) external potential is that of Eq. ([T]) 
where fi = — r is the dipole operator. 

Starting with an initial guess for the electric field e{t), 
we maximize the expectation value of the target operator 
6: 



Ji[^] = (*(r,r)|6|*(r,r) 



(6) 



where O = |$f) (*1'f| is now a projection operator, since 
we aim at maximizing the occupation of the target state 
<I>F at the end of the field at time T: 



Ji = |(vI/(r,r)|<i>F)| 



(7) 



In the following, this quantity is referred to the yield. 

As a constraint, avoiding fields with very high energy, 
the fluence (time-integrated intensity) of the field is lim- 
ited by a second functional. 



dte^(t) - Eq 



(8) 



where Eq is the fixed fluence [see Eq. (|13p below] and a 
is a time-independent Lagrange multiplier,^^ 

Finally, the satisfaction of the time-dependent 
Schrodinger equation [Eq. introduces yet another 
functional. 



J3[e,*,x] = -2Im / (xmdt^Himit)), (9) 

JO 

where x(t) is a time-dependent Lagrange multiplier. 

Variation of J = Ji -I- J2 -I- J3 with respect to x, 
and a lead to the control equations 



e{t) 



H{t)^{t), *(0) = (10) 
H{t)x{t), x{T) = d^iT), (U) 

(12) 



I. 



—lm{xitM^{t)) 
a 



dte^{t) = Eq. 



(13) 



III. METHOD 

In OCT the objective is to find an external time- 
dependent field e{t) that drives the system into the prede- 
fined state through the solution of the Schrodinger equa- 
tion, 



i^^^ir,t) = H[ekit)]^{r,t). 



(4) 



Here e(f) is an electric field (e.g., laser pulse) dealt with 
the dipole approximation, so that the Hamiltonian has 
the form. 



(5) 



which can be solved iterativelyi^i^^ We apply a numer- 
ically efficient forward-backward propagation scheme in- 
troduced by Werschnik and Gross.^^ When solving the 
control equations, the Lagrange multiplier a is calculated 
through the fixed fluence Eq as explained in detail in 
Ref. I21I . The field is constrained by an envelope function 
of a form 



fit) 







a 








|Erf 




H)] 


+ Erf 
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T 







T 

T y ' ' ~b 

(14) 

with a = 100 and b = 20. This corresponds to a 
step function ascending (descending) rapidly at i ^ 6/4 
{t ^ T ~ b/A). The scheme also allows straightforward 
inclusion of spectral constraints discussed in the follow- 
ing section. In the numerical calculations we have used 
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the octopus code'^*' which solves the control equations in 
real time on a real-space grid. 

To approximate the time-propagator we applied the 
time-reversal symmetry, i.e., propagating 5'(t) forward 
by Ai/2 should correspond to propagating 5'(t -f At) 
backward by At/2. This condition leads to an approx- 
imation for the propagator, ■^-'^ which can be further im- 
proved by extrapolating the time-dependent potentials. 
In octopus^ the used method is called Approximated 
Enforced Time- Reversal Symmetry (AETRS). 



IV. RESULTS 

First, the system is prepared from the ground state |0) 
to the desired superposition. Hence, we simply define 
the target wave function in Eq. (jH]) as |$f) = 1^)- Wc 
set the field length to T = 100 (~ 6 ps) and the initial 
frequency to 0.5 corresponding to the oscillator frequency 
LjQ of the DQD. Unless stated otherwise, the fluence is 
fixed to Eq = 0.3, so that the average intensities are of 
the order of 10"^ W/cm'^ (note the units given in Sec. HIl). 

The optimized field in Fig. [Tfa) looks rather compli- 
cated with distinct high-frequency components, whose 
role and possible removal is discussed in detail below. 
The occupations of the states, i.e., their overlaps with the 
time- propagated wave function, are plotted in Fig. [TJb) . 
The ground state (initially occupied) and the first ex- 
cited state (initially empty) get half populated, so that 



FIG. 3: (Color online) (a) Optimized fields without (thin blue 
line) and with spectral constraints (thick red line) for transi- 
tion |L) — >■ \R). (b-c) Occupations of the five lowest eigen- 
states during the process, (d) Occupations of the initial and 
target superposition states. 



their superposition \L) becomes fully populated and the 
electron is localized in the left well. The obtained yield 
is as high as 0.99985. 

After the preparation of the localized state we optimize 
a transition from \L) to \R), i.e., a charge transfer be- 
tween the quantum wells. The result of the optimization 
is summarized in Fig. |31 The optimal field having a fixed 
duration of T = 100 [thin blue line in Fig. Ufa)] leads to 
an extremely high yield of 0.9992. In Fig. [Sfb) and (c) 
we plot the occupations of the five lowest states during 
the charge-transfer process. Each of these states reach a 
maximum occupancy of more than 10 % during the pro- 
cess. The tenth lowest state still obtains '-^ 1 % of the 
occupation. Thus, with the present length of the field, 
the inclusion of several states seems to be crucial for the 
success of the optimization. Consequently, an alternative 
OCT procedure for a few-level model system (higher lev- 
els omitted) would lead to a completely different solution 
field, which most likely would perform poorly when ap- 
plied to the "full" system (as here) due to the leaking of 
the occupancy to higher states^ 

Similarly to the preparation field in Fig. HJa), the op- 
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0.9999 



FIG. 4: (Color online) (a) Spectrum of the optimized field 
for the process \L) ^ \R) without (blue thin line) and with 
a spectral constraint (red thick line) at cuth = 0.817. (b) 
Occupation of target state as the function of the frequency 
threshold cjth used as the filter. The step-like form of the 
curve is a consequence of the discrete Fourier transform. 



timized charge-transfer field in Fig. IH^a) shows abrupt 
peaks corresponding to high frequencies. Hence, the field 
would be practically impossible to construct, e.g., with 
the present pulse-shaping techniques. To relieve these 
limitations, we apply a spectral constraint cutting off the 
high-frequency components beyond a selected threshold 
frequency cjth- The thick red line in Fig. [3l^a) shows the 
field obtained using wth = 0.817 (~ 14 THz) in the opti- 
mization. The Fourier spectra of both fields are shown in 
Fig. Ufa). Both fields have a peak at at a; = 0.5, which in 
fact corresponds to the oscillator frequency ujq in Eq. ([ij . 

It is interesting to note that despite the relatively 
strong frequency constraint at Wth = 0.817, leading to 
a considerable smaller search space for the optimization, 
the obtained yield is reduced only down to 0.9986. This 
is a significant result in view of the fact that the original 
field has a large fraction of high frequencies as shown in 
Fig. Hl^a). Nevertheless, using a frequency filter does not 
considerably reduce the importance of higher states in 
the optimization: in this particular case the fifth lowest 
state still gains a maximum occupancy of ~ 10 %. In any 
case, further tightening of the threshold to smaller val- 
ues leads to decrease in the overlap as demonstrated in 
Fig. Hfb). The dependency is nonmonotonic due to nu- 
merical variation (note the logarithmic scale) and has a 
step structure resulting from the discrete Fourier trans- 
form. Below wth 0.5 corresponding to the oscillator 
frequency the fidelity collapses from 96 % to 42 %. If the 
fluence of the field is increased from 0.3 to 1, the critical 
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FIG. 5: (Color online) (a) Obtained yield for the process 
\L) — > \R) as a function of the field length, (b) Yield for 
the same process as a function of the fluence for three field 
lengths. 



threshold remains at 0.5, at which the fidelity decreases 
from 92% to 74%. 

Besides the threshold frequency, the main constraints 
in the field to be optimized are the length and the fiu- 
ence. In Fig.[S{a) we show the yield, again for the process 
\L) ^ as a function of the field length for fixed fiu- 
ence values Eq — 0.3 and 1, respectively. Both cases show 
some saturation around T > 100 although, as expected, 
the smaller fluence allows longer fields with even higher 
fidelities. However, increasing the yield above 0.9999 is 
difficult in this fiuence range unless relatively long fields 
are required. Here, the chosen length T — 100 seems 
an appropriate compromise between T and the obtained 
yield. 

Figure [SJb) shows the yield as a function of the fluence 
for three flxed field lengths. We remind that the fluence 
is a time-integrated quantity [see Eq. (jl3l) ] so that the 
curves correspond to different distributions of the energy 
in the field. In all cases the yield first increases expo- 
nentially with the fiuence until a point of saturation is 
reached. When T = 200 the slight decrease in the overlap 
at fiuences above ~ 0.2 might be due to numerical con- 
straints: in that regime higher and higher states (with an 
increasing number of nodes) are required, and they have 
a finite accuracy on the numerical grid. 
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FIG. 6: (Color online) Optimized field (upper panel) for the 
five- fold charge-switch process \L) ^ \R) ^ |L) ... — >■ |_R) 
(lower panel) in a double quantum dot. Here we have used a 
threshold frequency of uth ~ 0.817 leading to the target state 
occupation of 99.46 % at the end of the total five-fold process. 



TABLE I: Final-state occupations after n-fold sequential 
charge-switch processes obtained when merging the optimized 
fields. They are compared with an estimate based on the 
power of the yield given by the original (not inverted) field 
(see text). 



Eo uith n : 



5 n = 10 



50 



100 



0.3 



oo 0.9992(4) 
power law 



0.9989 
0.9985 



0.9975 
0.9962 



0.9961 
0.9924 



0.9848 
0.9625 



0.9899 
0.9263 



0.3 0.817 0.9985(9) 0.9977 0.9946 0.9865 0.9562 0.9589 

power law 0.9972 0.9929 0.9859 0.9317 0.8681 

0.3 0.629 0.9954(7) 0.9896 0.9606 0.8874 0.7092 0.7200 

power law 0.9910 0.9775 0.9556 0.7968 0.6349 



1.0 



oo 0.9990(3) 
power law 



0.9984 
0.9981 



0.9964 
0.9952 



0.9919 
0.9903 



0.9952 
0.9526 



0.9863 
0.9074 



1.0 0.817 0.9955(4) 0.9909 0.9842 0.9837 0.9838 0.9856 
power law 0.9911 0.9779 0.9563 0.7999 0.6398 



Finally we consider sequential charge-transfer pro- 
cesses by merging optimized fields together. For the pro- 
cess \L) — >• — >■ |L) — >■ . . . we combine, in turns, the 
optimized field cl^t?. (see above) with its time-inversion 
corresponding to ^r^l- The combined field with a 
threshold frequency wth = 0.817 is visualized in the upper 
panel of Fig. [51 The lower panel shows the occupations 
of the states \R) and \L) by solid and dotted lines, re- 
spectively. The final yield after the five-fold process is 
99.46 %. 

A more complete view on the results of up to 100 se- 
quential processes is given in Table ID We consider fiu- 
ences Eq = 0.3 and 1 for a single process, respectively, 
and different threshold frequencies as well as the case 
without a filter, i.e., Wth oo. The total yield shown 
in the table can be expected to (roughly) follow a power 
law, Ji^tot = "^"singiej where n is the number of processes 
(charge transfers). In this respect, the fidelity for a sin- 
gle transfer is essential for the quality of the final result. 
Indeed, the computational result follows the trend of the 



power law, but we find also significant differences: most 
importantly, in all cases the computational result is bet- 
ter than the prediction of the power law. The most dra- 
matic discrepancy can be found in the last example with 
Eq = 1 and Wth = 0.817, where after 100 pulses the yield 
is still almost 99 %, whereas the power law predicts is 
only 64 %. The reason behind the robustness of the yield 
in a sequential process is in the identity of the frequency 
components between the original and inverted fields, so 
that the population "lost" in higher states is partially at- 
tained back in the inverse process. There is, however, no 
clear trend in Table U indicating which field parameters 
are particularly favorable for robust sequential processes. 
Construction of such population-preserving, yet well op- 
timized sequential fields is a subject of future work. 

We point out that a critical aspect in the feasibility 
of the present approach is the sensitivity to decoherence. 
Typical decoherence mechanisms in semiconductor quan- 
tum dots are the hyperfine effects and interactions with 
optical and acoustic phonons. Their interplay and sig- 
nificance are largely dependent on the external condi- 
tions in a particular device. Detailed assessment of these 
mechanisms is beyond the scope of this work. We only 
mention that typical decoherence times in semiconductor 
quantum dots have been measured to be relatively large, 
even up to the millisecond scale, '^^ which in fact has been 
one of the main motivations of utilizing quantum dots in 
solid-state quantum computingi^i In view of the time 
scales considered here (up to hundreds of picoseconds) 
we believe that our approach is robust against the essen- 
tial sources of decoherence, although further analysis is 
in order. 



V. SUMMARY 

Here we have numerically constructed optimal fields for 
charge-transfer processes in single-electron double quan- 
tum dots. The only approximation has been the model 
potential for the device, so that no truncation of eigen- 
states in terms of A'^-level approximations have been used. 
We have found that optimal control theory provides an 
efficient way to operate on desired superpositions of the 
eigenstates regarding both the preparation of the local- 
ized state as well as coherent charge transfer between 
the quantum wells. We have analyzed the interplay be- 
tween different field constraints including the frequency 
filter, fluence, and the field length. Relatively strict fre- 
quency filters can be used without losing the extremely 
high yields obtained in the processes. Combination of 
the optimized pulses can be used in sequential charge 
transfers while preserving the high fidelity. 
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